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The large-distance asymptotic behavior of the field-field correlators has been computed for one- 
dimensional impenetrable anyons at finite temperatures. The asymptotic behavior agrees with the 
predictions of conformal field theory at low temperatures and reproduces the known results for 
impenetrable bosons and free fermions in appropriate limits. We have also obtained an integrable 
system of partial nonlinear differential equations which completely characterizes the 2-point corre- 
lation functions. The system is the same as for bosons but with different initial conditions. 
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Anyonic exchange statistics is defined naturally for 
impenetrable hard-core particles in the two-dimensional 
space [J, Q , which is sufficiently "large" to allow particles 
to rotate around each other, but, on the other hand, is 
"small" enough for the particle trajectories to be entan- 
gled by such rotations, so that the winding number of 
the rotations is a well-defined characteristics of each tra- 
jectory. By contrast, one-dimensional (ID) space seems 
too small to support any exchange statistics, since im- 
penetrable particles can not be exchanged at all when 
confined to move on a line. This simple expectation, fre- 
quently expressed in the literature, is not fully correct 
because of two circumstances. One is that in the case of 
quasi-periodic boundary conditions (i.e., motion on a cir- 
cle), the particles can be exchanged indirectly by trans- 
fer along the circle. These indirect exchanges mean that 
the statistics affects the precise form of the boundary 
conditions and, therefore, the position of the quasipar- 
ticle energy levels, making the thermodynamic proper- 
ties of the system of N impenetrable particles statistics- 
dependent at the level of 1/N corrections [!, 0, H, 
Even more importantly, in all realistic ID systems, par- 
ticles can be added to or removed from the system, pro- 
cesses which are described by the field-field correlation 
functions, and are strongly sensitive to the statistics even 
in the limit N — ► oo. Indeed, due to fluctuations in the 
ID quantum liquid, its multi-particle wavefunctions con- 
tain terms with different distribution of particles relative 
to a fixed point x, where the external particle is added. 
This means that the amplitudes of the addition/removal 
processes, and associated field-field correlators unavoid- 
ably involve contributions with different ordering of par- 
ticles, and therefore depend on the exchange statistics 
• The purpose of this work is to present exact calcu- 
lation of the long-distance asymptotics of the field-field 
correlators for the system of ID impenetrable hard-core 
anyons at arbitrary finite temperatures, < T < oo. In 
practice, the anyonic liquid studied in this work is simi- 
lar to the systems that can be realized with edges of the 
electron liquids in the FQHE regime. 



We consider the most basic model of ID anyons Q , the 
Lieb-Liniger gas, characterized by the second-quantized 
Hamiltonian 
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-htf(x)*{xj\, (1) 

where c > is the coupling constant and h is the chemical 
potential. We limit ourselves to the case of impenetrable 
particles, c — > oo. Anyonic exchange statistics is intro- 
duced by imposing the following commutation relations 
on the field operators: 



y{xi)^(x 2 ) 



^(x 2 )^(xi) 



(2) 



where k G [0, 1] is the statistics parameter, and e(x) — 
x/\x\, e(0) = 0. The commutation relations Q become 
bosonic for k = 0, and fermionic for k = 1. They imply 
that the multi-particle wavefunction acquires the phase 
factor e ±47rK upon the exchange of any of the two nearest- 
neighbor particles, where the signs ± need to be assigned 
consistently for all pairs of the wavefunction coordinates 
0, 0|- Thermodynamic properties of the Lieb-Liniger 
gas of anyons were studied in [!, 0, H, H, and for 
c — > 00, coincide (neglecting 1/N terms) with those of 
the free fermions. At T = 0, the quasiparticle momenta 
fill out the Fermi sea with the Fermi momentum and 
velocity given by, respectively, kp = Vh and Dp = 1\fh, 
in conventions of Eq. ([T]). Here, we are interested in the 
equal-time field-field correlator defined as usual by 



(3) 



The correlator depends only on the difference x = 
(xi — x 2 ) of the coordinates, and since (x 2 ) 1 $(xi)} — 
(xi)^(x 2 )}* , we can consider only the range x > 0. 
As the reference point for the results of this work, 
we mention that at T — > 0, the long-distance asymp- 
totics of the correlator (J3j> can be obtained from the 
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standard bosonization approach [10]. The particle op- 
erators are expressed in it through two bosonic fields: 
the integral 9(x) of the long-wave part of the density 
fluctuations p(x), i.e., d9(x)/dx = p(x)/w, and the con- 
jugate field 4>(x) defined by the commutation relation 
\4>{x),0(x')] — iire(x — x')/2. The usual bosonization 
ideology, i.e., the operator e 1 ^ reducing density by the 
amount that corresponds to one particle, and 9 produc- 
ing the appropriate phase changes of the wavefunction 
across each particle, implies immediately that the any- 
onic field operators ^ can be written as [ll| 

t) ~ e -i(K+2m)[k F x+8(x,t)] £ i<t>(x,t) _ ^ 

m 

The sign ~ in Eq. Q is the reminder that the relative 
amplitudes of different components are not defined in this 
sum. The commutation relations between 9 and <p mean 
that the individual terms in Q satisfy the necessary ex- 
change relations |2]) away from xi = X2- Using Eq. ((4]) 
for the field operators, and taking the standard average 
over the equilibrium fluctuations of 9 and 4>, one finds 
the long-distance behavior of the correlator ([3|): 
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[sinh( 7 rTx/^)]( 1 + K2 )/ 2 



-l e 
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where the constants b are some undetermined amplitudes. 
Equation (O gives the two leading terms which corre- 
spond to m — 0, — 1 in (|3} and are relevant for the dis- 
cussion below. In general, the correlator contains the 
hi gher -order terms which include the full expansion (j4]) 



111 ] and goes beyond it, if one uses the somewhat more 



general conformal field-theory approach Q 

The main new result obtained in this work is the ex- 
pression for the asymptotics of the anyonic field-field 
correlator ((3]) for large distances, x — ► +oo, that ex- 
tends the bosonization result to arbitrary temperatures. 
The bosonization ([4]) is known to rely effectively on ap- 
proximation of the linear quasiparticle spectrum: = 
±vp(k =F kp) for quasiparticle momenta k ~ ±/cf, which 
implies the exact electron-hole symmetry. This means 
that from the physics perspective, such an extension is 
needed for all phenomena (e.g., thermo-electric effects) 
which depend on the asymmetry between electrons and 
holes. On the other hand, our result is also interesting as 
an example of exact calculation of a non-trivial quantity 
in statistical mechanics. It is: 
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where cq, c_i are some complex constants, (3 = h/T, and 
C((3, k) is defined as 



C(/3,k) = - 

7T 



-Hoc. 



In I 



e x2 ~P + 1 

sA 2 — /3 plTTt 



dX. 
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The branch of the complex logarithm in Eq. {7J is fixed by 
the conditions that there are no cuts intersecting the real 
axis of A, and that the logarithm tends to zero at A — > oo. 
The complex constants Ao and A_i in ((6]) depend on (3 
and K as Xj = (— l) J [/3 + in(n + 2j)] 1 / 2 , or, separating 
explicitly real and imaginary parts, 

A, = [(-l)'(a, + (3) 1 / 2 + i( aj - (3) 1/2 }/V2, (8) 

where (Xj = \(3 2 + tt 2 (k + 2J) 2 ] 1 / 2 . Strictly speaking, the 
accuracy of our calculation allows us to keep the second 
term in ^ only when it does not decay much faster than 
the first: 9A_i < 2SAo, condition satisfied when the 
statistics parameter k is not too small. For n approach- 
ing the fermionic point n = 1, the two terms have the 
same decay rate 3Ao ~ S?A_i, and 5RAo ~ — 3?A_ i - see 
Eq. ©, so that the asymptotics exhibits beats with an 
exponentially decaying envelope. 

The method we used to obtain Eq. ([6]) (described be- 
low) allows us to find both the exponents and the am- 
plitudes Cj in this equation. More technical results for 
the amplitudes will be presented in the follow-up publi- 
cation. Here we focus on the main exponential factors. 
We start the discussion with the analysis of C(f3, k) ([7]) 
in different limits. In the "gas phase" , when the chemical 
potential h of the anyons is negative, i.e. j3 = h/T < 0, 
we can split the logarithm in Eq. (J7]) in two and use for 
them the expansion ln(l + z) — J2^=i(~ l) n+1 z n /n at 
\z\ < 1, to obtain: 



C(J3,k) = -= 



1 ^ 1 



'n ' — ' n 

n— 1 



3/2 I 
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i.e., for (3 — > — oo: 

C((3, K ) = (l/V^)e~ ll3l (l + e*™). 



(9) 



(10) 



When the chemical potential and (3 are positive, we di- 
vide the integration range in {7} so that for each interval 
one can use the same expansions of the logarithms, and 
find 



COM) = -\i-K{l-K)^+y j -([e 
tt L * — ' n \ 



-(-l)"]x 



71=1 



v^7 



z^'-fldX + [e m ™ - (-1)"] / e'^'^dX 

For large (3, both integrals in this expression are eval- 
uated asymptotically as \/(2nyf]3) + 0(/3 -3 / 2 ), and the 
resulting sums over n can be taken explicitly using the 
(6) standard formulae J2T=i (-1)" +1 /™ 2 = tt 2 / 12 , and 
Efc=i cos nTTn/n 2 — tt 2 B2(k/2), where B2(x) = x 2 — 
x + 1/6 is the second Bernoulli polynomial. This gives 
C(P,k) for (3 — > oo, i.e., in the low-temperature limit 
T < h, as 

C(/3, K ) = ^=(l- K ) 2 + 2z^(l-«). (11) 
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The point j3 = is regular in Eq. ([7]), and the expansion 
for small positive [3 — > +0 coincides with the expansion 
© for negative (3, when j3 — ► — 0. 

The regime of positive chemical potential and low tem- 
peratures, P — > oo, should correspond to the conformal 
limit, when the asymptotics of the correlator reduces to 
the bosonization result ([5]) . Using Eq. (TTTj) for the large- 
(3 behavior of C((3, k) together with Eq. © for Xj which 
gives in this limit: 



Ao = y/P + i 



TTK . i-x ,7r(2 — k) 



, (12) 



we see that the exact asymptotic behavior of the field- 
field anyonic correlator ([5]) reduces at low temperatures 
to the following expression: 



(*t(a;)#(0)) ~ c e -*(*T/«,)(« a +l)/2 e fa *JPK 



(13) 



•C_l e - a: ( 7r ' r / t 'F)[(«- 2 ) 2 + 1 ]/ 2 e lfc F(«- 2 ). 



The asymptotics (|5|) of the correlator within the 
bosonization approach agrees at non-vanishing temper- 
atures T / with the the exact expression (|T3|) . In- 
deed, for T 7^ 0, only one, the leading, exponent in 
saih^T x / v f) dominates at sufficiently large x, and the 
two expressions coincide. This means that, taken to- 
gether, the bosonization result and our Eq. (|13p de- 
scribe the correlator for all finite temperatures starting 
with T = 0, when Eq. ([5]) reduces to a combination of 
powers of x. 

Another important case reproduced by Eq. ([6]) is the 
bosonic limit k — > 0. As discussed above, only the first 
term in ([6]) needs to be taken into account in this case. 
The asymptotics depend essentially on the sign of the 
chemical potential h, i.e. (3. For (3 < 0, the argument of 
the logarithm in Eq. is positive and 



1 f + °° 

C(p, k = 0) = C{8) = - I In 

7T 



1 



o\*-f> _ 1 



dX. 



For /3 > 0, the argument of the logarithm is negative in 
the interval A € [— y/J3~, \ffi\, so that 



lim C(/3,k) =C(/3) + 2iJj3. 

K — >0 



Also, at K = Eq. §8§ gives Ao = \/~j3, and A = iy\P\ for 
positive and negative (3, respectively. Combining these 
equations, we see that Eq. ([6]) reduces at n — ► to 



*P^Cf(/ t /T)+ v '|/ l |] ft < q 



(*t(a,)*(o)) ~ c 



h>0, 



e -x^C(h/T) 



(14) 



and agrees with the results obtained previously for the 
impenetrable bosons (23. 24 1. 

In the opposite fermionic limit k — > 1 , all characteris- 
tics of the system (fTJ), including the correlators, should 



be those of the free fermions. Indeed, Eqs. ||7J and (JS]) 
show that C(/3, k) = and A = -(A_i)* = (j3 + in) 1 / 2 
for k == 1, and Eq. ([5]) reproduces the asymptotics of 
the field correlator for free fermions. More generally, us- 
ing the standard Wigner-Jordan transformation, one can 
express the anyonic fields @ of arbitrary statistics n in 
terms of (in our case free) fermions £ as 

*(ar,*) = e i < 1 ~ K ^ x ' t ^(x,t), n{x,t) = [ dx'p(x',t), 



where p = is the particle density operator. This rep- 
resentation provides a simple and useful interpretation of 
the C{j3, k) part of the asymptotics ([6]). In terms of the 
free fermions £, the anyonic correlator is: 

(*t(ar)*(Q)) = (e t (x)e" (K - 1) ™^(0)), (15) 

where n& = n(x) — n(0) is the operator of the number 
of particles in the interval [0, x]. This equation shows 
that the decay of the correlator is caused in part by 
the fluctuations of the number of particles rid creating 
the fluctuations of the accumulated statistical phase of 
the correlator. The average taken over these fluctuations 
leads to the correlator suppression. For large distances 
x, one can treat the fluctuations of simply as fluc- 
tuations of the number of free fermions in the system 
of size x. They are governed by the Fermi distribution 
r)(k) = l/le^-^/T + l], which gives the probability #(fc) 
for the particle with momentum k to be present in the 
system. Depending on whether each particle is in the 
interval [0, a;] or not, its contribution to the phase factor 
e m(K~i)n d j n correlator (fT5j) is either e^" -1 ' or 1. 
With this understanding, 

( e «r(*-l)n d ) = JJ[! _ + e W"-l)^(jfe)] 



exp 



f- 

l2vr 



dk In ■ 



,(k 2 -h)/T _ pi-KK 



,(k 2 -h)/T + 



r }- (16) 



The change of variables from k to A shows that this is 
exactly the C((3, k) part of the asymptotics ©■ 

Finally, we present an outline of the calculations lead- 
ing to the asymptotics ([6]) of the anyonic field-field cor- 
relator. The starting point of these calculations is the 
representation of the correlator as the Fredholm deter- 
minant of an integral operator which we have obtained 
previously by two different methods: the anyonic gen- 
eralization 14, [l5| of Lenard's formula 13], and direct 
summation of the form- factors [l6j]. The kernel of the 
integral operator of interest can be written as 



K T (X,p) 



e+(A)e_(,u) - e_(A)e + (M) 



2i(X - p) 



(17) 



where e±(A) = y/rJ(X) 



Here, and in the remainder 



of the paper, to simplify the notations we use the co- 
ordinates rescaled by the temperature-dependent factor: 
x = (x\ — x 2 )Vt/2. 
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The integral operators with the kernels of the form 
([TTf represent a special type of the general class of "inte- 
grable" operators, which characterize a large set of corre- 
lation functions of integrable models and spectral distri- 
butions for random matrix ensembles 17, 1 



22j. An important property of the integrable operators 
is that the resolvent of the operator is also integrable. In 
the case of the kernel ([T7|) . the resolvent kernel can be can 
be written in the same form through the two functions 
f± (A) introduced as solutions of the integral equations 



/±(A)-7 



+ 00 

-oo 



Xr(A,/i)/±0*)d/i = e±(A) ) (18) 



where 7 = (1 + e 47rK )/7r. In terms of these functions, one 
can define the auxiliary "potentials" Bi m as 



B lm {x,(3,n) = 7 



e;(A)/ ro (A) dX, l,m = ±. (19) 



With these notations, the result for the determinant rep- 
resentation of the field correlator derived in 15, [l6| is 
expressed as 

(¥( Xl )*(x 2 )) =B ++ (x,P, K )e^ x -^ , (20) 

where a = lndet(l — -fKr) is the logarithm of the Fred- 
holm determinant. 

Followin g t he method developed for impenetrable 
bosons 21, 23, 24]], we have obtained the integrable sys- 
tem of partial nonlinear differential equations for the po- 
tentials B ++ and B^ : 



d?B + _ 



ld x dpB + ^ 



dxBj — — B , 



(21) 



with the initial conditions at fixed (3 and x — > 0: 



B ++ (x, (3, K ) = B+_{x, (3, ft) = jd(P) + [jd(l3)] 2 x , (22) 
where d{(3) = J 00 $(A)dA, and the condition 

B++(x,(3, k) = -B_| (x, (3, ft) = , for f3 — > — 00 . 

The derivatives of the logarithm of the determinant a are 



cter = -B, 



dia = -Bi 



dpa = -xdpB + _ + l -(d fi B + ^) 2 - l -{d fj B ++ f,(2?,) 

which shows that the field-field correlator can be ex- 
pressed in terms of the solutions of the integrable system 
([2"Tj) . Also, a satisfies the equation 

{dpd 2 x a) 2 + A{d 2 x (T)[2xdpd x <j+{di3d x a) 2 -2dp(j] =0 (24) 

with initial conditions 



o = ->yd(0)x - [ 7 d(/3)] 2 
a = , j3 — > —00 . 



0(x 3 







The integrable system (|2"Tj) is the same system that char- 
acterizes impenetrable bosons, but with different initial 
conditions (12"2"1) . At ft = 0, i.e. 7 = 2/7T, it reproduces 



the bosonic result [2l|, |23j, |24|. The same phenomenon 
was observed by Santachiara and Calabrese [25j (see also 
(26|) who showed that the field correlator at zero tem- 
perature in the finite box satisfies the same Painleve VI 
differential equation characterizing impenetrable bosons 
but with statistics dependent initial conditions. 

As the final step, the large distance asymptotics of 

the potentials -B++ , B^ were found from the matrix 

Riemann-Hilbert problem associated with the integrable 
system ([2~T]) . The technique is similar to the one devel- 
oped for impenetrable bosons 21, HH. The asymp- 
totic behavior of the potentials, combined with Eqs. ([2TJ]) 
and ([23]) . gives the main result of this work, the large 
distance asymptotics ((6|) of the field-field correlator for 
impenetrable one-dimensional anyons. 
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